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I. INTRODUCTION 

The study of the phase diagram of QCD is an im¬ 
portant theoretical challenge with numerous phenomeno¬ 
logical implications, e.g., for heavy-ion collision experi¬ 
ments, early Universe cosmology, and astrophysics. The 
properties of the theory along the temperature axis have 
been intensively studied by means of lattice calculations 
and it is by now well established that the first-order 
confinement-deconfinement phase transition of the pure 
gauge SU(3) theory becomes a crossover in the presence 
of light dynamical quarks with realistic masses |l| . 

The situation is much less under control at finite bary- 
onic chemical potential, where lattice simulations suffer 
from a severe sign problem [2, Q. Calculations with re¬ 
alistic quark masses are limited to the region of small 
chemical potentials in units of the temperature. One 
hotly debated issue in this context is the possible exis¬ 
tence of a line of first-order chiral transition with a criti¬ 
cal endpoint, as predicted by various model calculations 
0-[6:]. Efforts to tackle this question directly from the 
QCD action combine the investigation of various ways to 
circumvent the sign problem on the lattice 000 and 
the use of nonperturbative continuum approaches 0, 
where the sign problem, although not completely absent, 
is much less severe. 

Unlike lattice calculations, continuum approaches rely 
on identifying efficient approximations for the relevant 
dynamics. Standard perturbative tools are notoriously 
inadequate in this context because the typical transition 
temperatures are of the order of the intrinsic scale of the 
theory, where the coupling becomes large. This motivates 


the use of nonperturbative approaches, typically based 
on truncations of the functional renormalization group or 
Dyson-Schwinger equations fl0l - fl6l |. Lattice simulations 
in situations where they are well under control provide 
then valuable benchmarks for testing the various meth¬ 
ods and approximations. 

Examples of such situations are the case of a purely 
imaginary chemical potential, where there is no sign 
problem, or the study of QCD with heavy quarks. A sys¬ 
tematic expansion around the infinite mass (quenched) 
limit can be devised, for which the sign problem is mild 
enough so that numerical simulations can be performed 
up to large chemical potentials p~7l420(| . Although these 
do not correspond to physically relevant situations, the 
phase diagram in these cases is interesting in its own 
and actually features a rich structure. This has been in¬ 
vestigated both on the lattice and by means of effective 
models [2ll428j . 

In the present work, we investigate the phase diagram 
of QCD with heavy quarks by means of a modified per¬ 
turbative expansion in the context of background field 
methods. This is based on a simple massive extension of 
the standard Faddeev-Popov (FP) action in the Landau- 
DeWitt gauge which has been successfully applied to the 
confinement-deconfinement transition of the SU(2) and 
SU(3) Yang-Mills theories |29M3ll | . 

The motivations for such a massive extension are 
twofold. First, gauge-fixed lattice calculations of the 
vacuum Yang-Mills correlators in the (minimal) Landau 
gauge have shown that the gluon propagator behaves as 
that of a massive field at small momentum, whereas the 
ghost remains massless. This suggests that the dominant 
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aspect of the nontrivial gluon dynamics (in the Landau 
gauge) might be efficiently captured by a simple mass 
term, up to corrections that can be computed perturba- 
tively. This scenario has been put to test in Refs. 1321. ~ll . 
where a one-loop calculation of the gluon and ghost prop¬ 
agators in a simple massive extension of the Landau 
gauge—known as the Curci-Ferrari model [35) —has been 
indeed shown to describe the lattice results remarkably 
well. Similar successful results have been obtained for 
the three-point correlation functions [36| as well as for 
the vacuum propagators of QCD [37] and the Yang-Mills 
propagators at finite temperature [38| • An interesting as¬ 
pect of this approach is that the gluon mass acts as an 
infrared regulator and perturbation theory is well defined 
down to deep infrared momentum scales [341], 

The second motivation is more formal. It is related to 
the fact that the FP quantization procedure, which un¬ 
derlies standard perturbative tools, is plagued by the ex¬ 
istence of Gribov ambiguities [39| and is, at best, a valid 
description at high energies. For instance, it is known 
that a nonperturbative implementation of the BRST 
symmetry of the FP action on the lattice leads to unde¬ 
fined zero over zero ratios for gauge-invariant observables 
[13, EH- In fact, existing gauge-fixing procedures on the 
lattice—e.g., the minimal Landau gauge—typically break 
the BRST symmetry. This suggests that a consistent 
quantization procedure, which correctly deals with the 
Gribov issue, is likely to induce effective BRST breaking 
terms. The simplest such term consistent with locality 
and renormalizability is a gluon mass term . 1 

In Ref. [29[, we have extended this approach to the 
Landau-DeWitt gauge—which generalizes the Landau 
gauge in the presence of a nontrivial background field 
with the aim of studying the confinement-deconfinement 
phase transition of SU(fV) Yang-Mills theories. Back¬ 
ground field methods allow one to efficiently take into 
account the nontrivial order parameter of the transi¬ 
tion [43|. We have shown that a simple one-loop calcula¬ 
tion correctly describes a confined phase at low tempera¬ 
ture and a phase transition of second order for the SU(2) 
theory and of first-order for the SU(3) theory, with tran¬ 
sition temperatures in qualitative agreement with known 
values. Two-loop corrections, computed in Refs. 30l[3l[ . 
make this agreement more quantitative; see Ref. |44( for 
a short review. We mention that two-loop contributions 
also resolve some spurious unphysical behaviors (e.g., 
negative entropy) observed in ba ckg round field calcula¬ 
tions at one-loop order; see, e.g., [33j . 

The present work is a natural extension of these stud¬ 
ies to QCD with Nf = 2 + 1 heavy quarks flavors at 


1 A more precise relation between the Gribov problem and the 
gluon mass term has been obtained in Ref. 1421 . where a new 
quantization procedure based on a particular average over the 
Gribov copies along each gauge orbit has been proposed. The 
bare gluon mass is related to the gauge-fixing parameter which 
lifts the degeneracy between the copies. 


finite chemical potential. We compute the background 
field potential at one-loop order, from which we can read 
the value of the order parameters—the averages of the 
traced Polyakov loop and of its Hermitic conjugate—as 
functions of the temperature and of the chemical poten¬ 
tial. We show that this simple calculation accurately 
captures the rich structure of the phase diagram, both 
for real and imaginary chemical potential. Moreover, we 
obtain parameter-free values for the dimensionless ratios 
of the quark mass over the temperature at various critical 
points which compare well with lattice results. 

The paper is organized as follows. We present the basic 
framework, i.e., the formulation of the QCD action in the 
massive extension of the Landau-DeWitt gauge in Sec. HD 
and discuss in detail various symmetry properties of the 
relevant generating function in Sec. Hill This generalizes 
known material to the case of a nontrivial background 
field. The one-loop calculation of the background field 
potential is straightforward and is detailed in Sec. [IV] 
The rest of the paper is devoted to the discussion of our 
results at vanishing (Sec. [V]), imaginary (Sec. IVD) . and 
real (Sec. IVIII) chemical potential. In the latter case, we 
discuss how the sign problem also affects continuum ap¬ 
proaches. Finally, Appendix 1X1 briefly presents some con¬ 
sequences of charge conjugation symmetry at vanishing 
chemical potential and Appendix [B] provides an approxi¬ 
mate calculation of the background field potential which 
allows for a simple analytic understanding of some results 
presented in the main text. 

II. THE QCD ACTION IN THE MASSIVE 
LANDAU-DEWITT GAUGE 

The Euclidean action of QCD in d dimensions with N 
colors and Nf quark flavors reads 

N f 

J F H» F llv + Y2' t i ; f(V + M f+m o)Vv}> (!) 

/=i J 

where f x = dr f d d ~ 1 x, with (3 the inverse temper¬ 
ature, F« v = d^Al - d v A\\ + gf abc A b ^A c v , with g the 
coupling constant and f abc the structure constants of 
the SU(iV) group, and V^ip = — igA^t a ) ip, with 

t a the SU(iV) generators in the fundamental represen¬ 
tation, normalized as tit a t b = 5 ab /2. Finally, g de¬ 
notes the chemical potential. We leave the Dirac and 
color indices of the quark fields implicit and ip and ip 
are understood in the common sense as column and line 
bispinors, respectively. The Euclidean Dirac matrices 2 
7 ^ are Hermitian and satisfy the anticommutation rela¬ 
tions { 7 ^, 7 ,,} = 25^. 


2 These are related to the standard Minkowski matrices as 70 = 
7^ and 7, = -++. In the following, we work in the Weyl basis, 
where 7* 2 = 7^ 2 = 70,2 and 7*3 = 7^3 = -71,3. 
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The gauge field is decomposed in a background field 
and a fluctuating contribution as 

a; = a; + o,;, ( 2 ) 


t 8 = A 8 /2, where A a are the Gell-Mann matrices. Ac¬ 
cordingly, we consider the generating function (we make 
explicit the dependence on the chemical potential for lat¬ 
ter use) 


and the Landau-DeWitt gauge is defined as 

(D^r = 0, (3) 

with D ab = S ab d M + g/ oc6 j4° the background covariant 
derivative in the adjoint representation. The correspond¬ 
ing Faddeev-Popov gauge-fixing action reads 

S FP = J {(^cnzvr + ih a {D^) a } , (4) 

where D ab = + gf acb A\c and c are anticommut¬ 

ing ghost fields, and h is a Lagrange multiplier. Finally, 
as discussed in the Introduction, we also consider a mas¬ 
sive extension of the Landau-DeWitt gauge, with mass 
term 3 [2f| 


S m = J ^m 2 a“a“ . (5) 

In the following, we grab together the fluctuating fields 
of the pure gauge sector, which all belong to the adjoint 
representation of the gauge group, as tp = (a M , c, c, h). 
The total gauge-fixed action Sqcd + <Sfp + S m is invari¬ 
ant under the combined transformation of the fluctuating 
fields 


Z(J,A,n)= f , ( 8 ) 

with the source term 

sWc = [ j;a; (9) 

J X 

restricted to a constant source in the temporal direction 
and in the Cartan subalgebra: J“(i) = 5 ll oJ a , with a = 
3, 8. In the following, J and A denote vectors in the 
plane spanned by the SU(3) Cartan directions (to which 
we shall refer as the Cartan plane in the following): 

M!)- (10) 

The main quantity of interest in the following 
is the Legendre transform —In Z {j(A, A, p), A, gf) + 
f x J k (A, A, g)A k , where J k {A,A,g) is defined through 


/3tt dJ k 

We shall evaluate this Legendre transform at A = A, that 
is, 



if -> ipU ] , c p a t a -> Ucp a t a U t (6) 

and of the background field 

->■ UA“t a U f + -Ud^, (7) 


where U is a local SU(iV) matrix. 

The background field A is nothing but a gauge-fixing 
parameter and is used here as a device to capture non¬ 
trivial physics in an efficient way. We constrain it so 
as to break as few symmetries as possible. The sym¬ 
metries of the Euclidean space with the boundary con¬ 
ditions implied by finite temperature field theory allow 
for a constant vector field in the temporal direction: 
A“(*) = <S ;l oA“. Moreover, it is always possible, through 
a global color transformation, to bring the color vector 
A a in the Cartan subalgebra of the gauge group, spanned 
by the diagonal generators. In the following, we denote 
the latter by t k . For SU(3), these are t 3 = A 3 * * /2 and 


v (r,») = ~]nZ(j(A,fj),A,n) + J k (A,^)A k , ( 12 ) 

where r k = A k /{fig), D is the spatial volume, and where 
we defined J k {A,g) = J k {A, A, g). A priori, one could 
leave A independent of A and extremize the effective po¬ 
tential with respect to A at fixed A. However, one obtains 
the same physics by identifying A = A and extremizing 4 
with respect to A. This is a particularly convenient choice 
since the fluctuating gauge field has vanishing expecta¬ 
tion value, (a) = 0, even at nonvanishing source. 

Finally, important physical observables to be consid¬ 
ered below are the averages of the traced Polyakov loop 
in the fundamental representation and of its Hermitic 
conjugate: 



(\P exp 



Jo dTAa ° ta ) 




(13) 

(14) 


3 In principle, one could implement different bare masses for the 
gluon fields in the temporal and spatial directions, e.g., to ac¬ 
count for Gribov ambiguities on asymmetric lattices at finite 
temperature. We chose here the minimal extension of the FP 
theory with a symmetric, vacuumlike mass term. 


4 We shall see below that, for any imaginary chemical potential, 

one needs to minimize the potential G 3 - However, for real chem¬ 

ical potential, some amendments need to be made to this rule. 
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where P and P denote path ordering and antiordering 
respectively. In general, £(p) is not equal to the complex 
conjugate £*(p), e.g., in the case of a complex action. 
For real chemical potential, the quantities (fl3l) and ED 
are real and are related to the free energies F q and F q 
of a static quark or antiquark as £(p) oc exp(— (3F q ) and 
£(p) oc exp(— (3F q ) [3, [45ij. The physical loops £(p) and 
£(p) are to be evaluated at vanishing sources, that is, 
at an extremum of the background field potential (fl2l) . 
In order to discuss symmetries below, it is useful to in¬ 
troduce the following functions of the background field, 
defined at nonvanishing source J(A 1 p) A 0, 



where it is understood that the averages on the right- 
hand side are evaluated at A = A, that is, at (a) = 0. The 
physical loops £(p) and £(p) are obtained by evaluating 
the functions m and ED at the relevant extremum of 
the potential (fl2l) . 


III. SYMMETRIES 

In this section, we discuss various transformation prop¬ 
erties of the generating function ([HD, background field po¬ 
tential ED and Polyakov loops functions (fl5l) and (flbl) . 
For later purposes, it is useful to consider the general case 
where J k , A k , and p are complex numbers. We first dis¬ 
cuss the effect of charge conjugation and of complex con¬ 
jugation which relate the cases /i -<A — /i and /i -<A p*, re¬ 
spectively. Next, we analyze the consequences of the con¬ 
tinuous global and local color transformations. Finally, 
we discuss the Roberge-Weiss symmetry which relates the 
theories with chemical potentials p and p + 2in/(3(3). 


A. Charge conjugation 

We perform the change of variables 

ip -a Cp , , p -A Cp (17) 

under the functional integral ©. Here the charge con¬ 
jugation matrix C only acts on the Dirac indices of fields 
in the fundamental representation whereas C acts on the 
color indices of adjoint fields. The former satisfies 

= -it ( 18 ) 

and is given by C = 7072 hr the Weyl basis. It is such 
that C~ l = = C 4 = — C. As for the latter, it acts on 

the generators of the color group as 

C ab t b = _ (. t ay ( 19 ) 


Using the Gell-Mann basis, one easily checks that it is 
a diagonal matrix with eigenvalues +1 on the lines 2, 5 
and 7 and —1 otherwise. It is thus an orthogonal 0(8) 
matrix such that O” 1 = C* = C and which conserves the 
structure constants (C aa 'C bb ' C cc ' f a ' b ' c ' = f abc ). Notice, 
however, that it has determinant —1 and it is, therefore, 
not a color transformation in the adjoint representation. 

It is an easy exercise to check that the only effect of 
the change of variables ED on the generating function 
© is to change 

J —y C J — — J , A — y C A — — A , and p —^ — p , 

(2°) 

where we used the fact that C kk = — 5 kk in the Cartan 
subalgebra. We conclude that 

Z(J,A,p) = Z(-J,-A,-p). (21) 

This implies the relation J(— A, — A, — p) = —J(A,A,p) 
for the source defined in Eq. ED- It follows that 
J(-A, -p) = —J(A,p) and thus 


V(r,p) = V(-r,-p). (22) 

Similarly, we deduce the relation 

£(r,p) =£(-r,-p). (23) 

B. Complex conjugation 

We now consider the change of variables 5 

ip —> K,ip , Pp -a , p^-Kp, (24) 

where the matrix K, = 7173 only acts on Dirac indices 
and is such that /C -1 = KA = /C 4 = —K and 

^1^ = it- ( 25 ) 

As for the matrix K acting on the color indices of the 
adjoint fields, it satisfies 

K ab t b = - ( t a )* . (26) 

Using the fact that the generators t a are Hermitian, we 
see that K = C. Upon the additional change of integra¬ 
tion variable h a -A —h a and the transformation 

J —> KJ* = - J* , A -A KA* = -A* , and p- > p*, 

(27) 

one recovers the original form of the generating function 
with all complex coefficients replaced by their complex 
conjugate, that is , 6 

Z(J,A,»)=Z*(-J*,-A*,p*). (28) 


5 This is similar to the /C-transformation discussed in Refs. [ 26 [[ 27 l| . 

6 Note that the integration variables are not affected by the com¬ 
plex conjugation. 
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Following similar steps as before, 
J(—A*,g *) = —J*(A,g) and thus 

this implies 

V(r,g) = V*(-r\g*). 

(29) 

One also shows that 


£{r,g)=t{-r\g *) 

(30) 


and similarly for £. 

C. Global color symmetry 

The gauge-fixed action S'qcd + Sfp + is invariant 
under global SU(3) transformations 

■0 — > 1Z4 1 1 4> —t , <p — > R<p, (3i) 

where 1Z is a global SU(3) matrix and R the correspond¬ 
ing SO ( 8 ) transformation 

R ab = 2tr (TZk a TZt b ) . (32) 

Of particular interest are those global color transfor¬ 
mations which do not mix the Cartan subalgebra with 
the other elements of the Lie algebra (i.e., R is block- 
diagonal), such that they leave the source J and the 
background A in the Cartan subalgebra. For instance, 
the SU(3) matrix 

(0 1 0 \ 

Tlx = -1 0 0 (33) 

V° 0 V 

induces such a block-diagonal color rotation R\ whose 
restriction to the Cartan subalgebra reads 

i?f art an =diag (_ 1;1 ). (34) 

In the Cartan plane, this corresponds to a reflexion about 
the 8 -axis. We can therefore write 

Z (J, A,/i) = Z (R\J, RiA, g') . (35) 

As before, this implies that the source J(A,g) trans¬ 
forms covariantly, i.e., J(R\A,g) = RxJ(A, g), from 
which we conclude that 

V(r,g) = V(R ir ,g). (36) 

Similar considerations yield, for the functions (fT51) and 

£(r, g) = £(Rxr, g) (37) 

and similarly for l. 

It is easy to check that the SU(3) matrices 

/0 0 -l\ /1 0 0 \ 

77 2 = 0 1 0 and Hi = 0 0 -1 (38) 

\10 0 / \0 1 0 j 


also generate color rotations in the adjoint representation 
such that the Cartan components do not mix with other 
color directions. In the Cartan plane, they correspond to 
mirror images about two axes passing through the origin 
and making an angle of ± 7 t /6 with the 3-axis. Together 
with Rx , they generate all the color rotations under which 
the Cartan subalgebra is stable. These form a C^ v group 
which also contains rotations by an angle ±27 t/ 3 around 
the origin. This is the Weyl group of the su(3) algebra. 


D. Background gauge invariance 

We now come to analyze the consequences of the invari¬ 
ance of the action under the local transformation ® 01- 
In order to maintain the background field constant, in 
the temporal direction, and in the Cartan subalgebra, we 
consider local color transformations of the form U(t) = 
e *O//3)0 t , w hich simply generate a translation of the 
background field in the Cartan plane: A —> A + (f>/(Pg). 
True gauge transformations—which leave physical ob¬ 
servables invariant—are periodic in time, U(/3) = 17(0), 
which implies that the eigenvalues of the matrix <j) k t k 
must be multiples of 2tt. This is solved by 

0=jlMl+j2M2, (39) 

with Ji,2 integers and 

= (40> 

Both the integration measure and the action at vanishing 
source in ® are left invariant by the transformation @- 
0 - The change of the source term is trivial and we get 

In Z(j,A + 4>/(Pg),n) = In Z (J, A, g) J k <j) k 

pn pn p g ' 1 ’ 

We deduce that J (A + cj)/(pg), g) = J(A,g) and thus 
that 

V(r,g) = V(r + (t>,g), (42) 

with 0 given in Eq. (1591) . One also shows that 

£(r, g) = £(r + (j>, g) (43) 

and similarly for £. 


E. Roberge-Weiss symmetry 

Following Ref. 0] , we now consider local SU(3) trans¬ 
formations of the type studied in the previous subsection 
with U (r) periodic in time up to a nontrivial element of 
the center Z 3 of the gauge group: U(P) = e ±2lv ^ 3 U(0). 
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This results in translations A —> A + j\e\ + j 2 e 2 of the 
background field with j\ 2 integers and 7 

ei= 2n (-1/V3) ’ e2=2n (-i/U) • (44) 

As is well known, the pure gauge action- including the 
FP and mass terms—is invariant under such twisted 
gauge transformations, whereas the source term trans¬ 
forms trivially, as in Eq. m- This symmetry is explic¬ 
itly broken by the antiperiodic boundary conditions of 
the quark fields in the fundamental representation. To 
cope for this, we perform the additional change of vari¬ 
ables 

ijj —> , ip — (45) 

All terms in the functional integral J8|) are invariant 
except for the quark kinetic term, whose variation, 
Tw V7oV : can, however, be compensated by a shift of 
the chemical potential in the imaginary direction. A sim¬ 
ilar analysis as in the previous subsections leads to the 
following identities for the background field potential 

{ 27 7 r \ 

V(r,n) = V[^r + e h n + —J . (46) 

where j = 1,2. As is well known, the traced Polyakov 
loop gets multiplied by an element of the center under 
the twisted gauge transformations considered here. This 
leads to 

e (r, /x) = e -2i7r / 3 £ (r + e,-, M + ^) , (47) 

Hr, /x) = (r + e„ M . (48) 

Observe that e\ — e 2 = u\ + u 2 , which shows that if 
we use the property (l46l) twice so as to add and sub¬ 
tract 2 z7t/( 3/3) to the chemical potential, we retrieve the 
original potential, up to a gauge transformation of the 
form (TE1) . Similarly, using the property (HTH) thrice so 
as to add three times 2 z7t/( 3/3) to the chemical poten¬ 
tial and using the relation 2ei + e 2 = u 2 we see that we 
generate a gauge transformation of the form (TT21) up to a 
translation by 2iir//3 of the chemical potential. We, thus, 
recover the well-known fact that the physics is unaffected 
by a change 8 jj, —> n + 2in/j3. 


7 The vectors ei and e 2 generate a lattice dual to that generated 
by the roots of the algebra; see below. This property generalizes 
to any compact Lie group with a simple Lie algebra, see [ 3 lll for 
a detailed discussion. 

8 This can be checked by performing the change of variables 

7/1 —y e 2 * 7 rr /^ , 0 , —y e~ 2 i 7 rr /^ 7 /;, which preserves the antiperi¬ 

odic boundary conditions. 


F. Vanishing sources 

We end this section by discussing the case of vanishing 
sources J = 0, which corresponds to the physical point. 
For instance, the partition function Z(fi ) is obtained as 
the generating function J5]) at J = 0. In the present 
context, where we set A = A, it is obtained from Eq. m 
at J = 0, namely, 

= (49) 

where r ex t(^) is an y physical extremum of V’(r,/x). In 
general, there are many such extrema, which form a set 
(Text(/•*)}■ The properties (1^1) and CTH) imply 9 that this 
set satisfies, 

{r-ext(^)} = {-7-ext (-/■*)} = (Vex t (/ 7 ± 2 x7t/(3/3)) =F ej} , 

(50) 

with j = 1,2. Similarly, the relation (TSUI) implies that the 
set of extrema of V * (r, /x*) is {-<xUsing Eq. (0D 
and the relations ([22| I . (l29l) . and (l46l) . we then retrieve 
the standard relations for the partition function Q 

m = Z(-V) = = z (/X + . (51) 

Similarly, using the definitions £(n) = l (r ex t (/x), /x) and 
= £(r ex t(/x),/x) and assuming that none of the 
symmetries mentioned here is spontaneously broken, we 
get the known relations for the averages of the traced 
Polyakov loop and of its Hermitic conjugate: 

m = U( M *) = e- 2 ™ l3 £ (m + ^) • (52) 

For imaginary chemical potential, the Roberge-Weiss 
symmetry can be spontaneously broken. In that case, 
the averaged Polyakov loops can be multivalued and the 
relations in (l5!?l) relate the sets of their possible values. 

IV. THE BACKGROUND FIELD POTENTIAL 
AT ONE-LOOP ORDER 

Here, we detail the calculation of the leading-order, 
one-loop effective potential (fl2l) . We consider a general 
compact Lie group with simple Lie algebra and we spec¬ 
ify to SU(3) when needed. For completeness we also give 
the corresponding leading-order (tree-level) expressions 
of the Polyakov loop functions (fl5l) and (usd. Our ap¬ 
proach is similar to that used in |23j . 


9 This is obviously true for the set of absolute minima which are the 
relevant extrema at vanishing and imaginary chemical potential. 
For a real chemical potential, the choice of the relevant extrema 
is more subtle, as discussed below. Our prescription, in Sec. [vm 
is also compatible with Eq. (Eg. 
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The one-loop contribution to the potential only in¬ 
volves the action up to quadratic order in the fluctuating 
fields. The pure gauge contribution has been computed 
in Ref. [29j. Denoting by M the quark mass matrix, the 
relevant contribution from the quark sector is 

V q {r, p) = —TrLn (ft + M + pj 0 - ig^ 0 A k t k ) , (53) 


simple Lie algebra. After taking into account the con¬ 
tribution from the massive gluon modes and the partial 
cancelation between massless modes contributions from 
the ghost and the h — a sectors, one has, in d = 4, 

^gauge(f) = 2'^' m ( r ) ~ (59) 


where the functional trace Tr involves a sum over 
fermionic Matsubara frequencies, an integral over spa¬ 
tial momenta and a trace over Dirac, color, and fla¬ 
vor indices. The Cartan generators t k can be diagonal¬ 
ized simultaneously and their respective eigenvalues p k 
form a set of dp (possibly degenerate, i.e., identical) vec¬ 
tors p in the space spanned by the Cartan directions, 
with the dimension of the fundamental representa¬ 
tion . 10 In group theory language, these are called the 
weights of the representation. The fundamental repre¬ 
sentation of the su(3) algebra has three nondegenerate 
weights, p € {pi, p 2 , P 3 }, which are related to the vec¬ 
tors ei j2 introduced in Eq. flU) as pi , 2 = —ei, 2 /( 47 r) and 
Pi + P 2 + P 3 = 0. Explicitly, 

Pl = 2 ( 1/V3 ) ’ 

P3 = ^(l/ 0 v^) ' 

Having diagonalized the color structure in Eq. (l53l) . we 
see that the constant temporal background A can be ab¬ 
sorbed in a redefinition of the chemical potential p —> 
p — ir p /ft , with r p = r k p k , which lifts the degeneracy 
between the various color states. We thus get 

V q (r,p) = Y V f(v ~ ir p/P)> ( 57 ) 

f,P 


(54) 

(55) 

(56) 


where the sum runs over all flavors and all color states 
(weights p) in the fundamental representation and Vjr ( p ) 
is the one-loop contribution from a single quark flavor 
/ in a definite color state at vanishing background field 
and nonzero chemical potential p [241 ] . This is given by 
the standard formula (written here for d = 4 and letting 
aside a trivial T- and /x-independent piece) 


V f °(p) 



1 + e ~d(e f q +M) 


(M — M)}j 


(58) 


where q 2 + Mj. 

The pure gauge contribution has been computed in [29] 
for SU(2) and SU(3). It is instructive to derive the cor¬ 
responding formula for a general compact Lie group with 


10 The present calculation easily generalizes to fields in any repre- 
sent at ion, provided one works with the corresponding weights. 


where 


J- m (r) = TrLn (D 2 + m 2 ) (60) 

and D p = d p — igS p oA k T k , with T k the generators of 
the Lie algebra in the adjoint representation. Here, the 
functional trace Tr involves a sum over bosonic Matsub¬ 
ara frequencies, a spatial integral and a trace over color 
indices. 

As for the quark contribution, this trace can be easily 
evaluated using the weights of the adjoint representation. 
In the subspace spanned by the Cartan directions, these 
form a set of dA vectors k whose components are the 
eigenvalues of the Cartan generators T k , where dA is the 
dimension of the adjoint representation. One sees that 
the temporal background enters as a shift of the Matsub¬ 
ara frequencies, which can be interpreted as an effective 
imaginary chemical potential —*r K //3, with r K = r k n k . 
We thus get 


rp r- OO 

dqg 2 In [1 - e-^ +ir “] , (61) 

where the sum runs over all color states (weights k) in 
the adjoint representation. By definition, the generators 
T k have as many zero eigenvalues -corresponding to k = 
0—as there are Cartan directions. Furthermore, one can 
show that the nonzero weights always come in pairs k = 
±a, which are called the roots of the Lie algebra. We 
obtain, 


T r°° r 

^ra(r) = ^2 J d( iq 2 \dc In [1 - e -/3e «] 

+ Yj In [l + e~ 2/3eq - 2 e~P eq cos r a ] j, (62) 

a 


where dc is the dimension of the Cartan subalgebra and 
the sum runs over the (dA — dc )/2 possible pairs (a, —a). 
For SU(3), there are two Cartan directions and dA = 8 so 
there are six roots a € {iaq,±c* 2 ,± 0 : 3 }, where we can 
choose [see Eq. m] 0 : 1,2 = «i,2/(47t) and 01+02 + 03 = 
0. Explicitly, 


a2 = K-V3)’ 

03 = -(j). 


(63) 

(64) 

(65) 






The sum in (l62l) runs only over the set {on, c* 2 , 0 : 3 }. For 
completeness, we mention that, for r a G [0, 27t], the func¬ 
tion jFm—o (r) admits the closed form [471 14§| 



O' 


Finally, the leading-order, tree-level expressions of the 
Polyakov loop functions m and (E2) are trivially ob¬ 
tained as 


£(r,n) = £(-r,p) 


tr (e irktk ) 
cLf 



(67) 


and do not depend explicitly on p. For SU(3), we have 


= - 


e *7% + cos(r 3 / 2 ) 


( 68 ) 


It is easy to check that the above expressions reproduce 
the SU(3) results of Ref. [2Ej. Also, it is a simple exercise 
to verify that the leading-order expressions m, 
m, and (1681) satisfy all the symmetry properties 

derived in the previous section. To this aim, it is useful 
to note that the scalar products m ■ pj , Ui ■ aj , and ej • aj 
are all multiples of 27r, whereas ei • aj = 2n/3 mod 2n 
and e 2 • aj = —2tt/3 mod 2n. 


V. VANISHING CHEMICAL POTENTIAL 

We now apply the above calculations to various situa¬ 
tions of interest. We begin our discussion with the case 
p = 0, where the lattice simulations are well under con¬ 
trol. We first discuss how the symmetry considerations 
of Sec. IIIII constrain the background field potential at 
p = 0. Then we study the phase structure of the theory 
as a function of the quark masses and we compare our 
findings with lattice results. 



FIG. 1 : The diagram summarizes the symmetries of the 
background field potential at vanishing chemical potential in 
the plane (r3,rg), with r 33 G R. The plain lines are mirror 
symmetry axes which reflect the global color symmetry and 
the background field gauge invariance of the original action, 
see Eqs. © and Q. The latter induces the translation sym¬ 
metries along the vectors ui,2, defined in Eq. (TIUll . These are 
symmetries of the potential also at nonvanishing p. For p = 0 , 
charge conjugation invariance results in an inversion symme¬ 
try about the origin and the dashed lines thus also become 
mirror symmetry axes. The white dots are all equivalent by 
combination of these symmetries. The shaded triangle repre¬ 
sents an elementary cell that is repeated all over the plane. 
In the pure gauge theory, the potential is also invariant un¬ 
der translations along the vectors ei,2, defined in Eq. (02), 
such that the black and the white dots are all equivalent by 
symmetry. The elementary cell is reduced to a third of the 
shaded triangle, i.e., to any of the three subtriangles shown 
with dotted lines. 


A. Symmetries 

A direct consequence of the relations (12T1) and (E51) at 
p = 0 is that the generating function is real if we choose 
the source and the background field components to be 
both real. Indeed, these relations imply 

Z(J,A, 0) = Z{-J,-A, 0) = Z*{J,A, 0). (69) 

In particular, this guarantees that the average value A 
of the gauge field at nonvanishing source, defined in 
Eq. (1111) . is real. This is important in the present setup 
where we eventually choose A such that A = A. More¬ 
over, the above relations, together with Eqs. (l23l) and 
(1301) imply, for a real background field, 

(70) 


and 

£(r, 0) = £(—r, 0) = £*(—r, 0). (71) 

In this case, the physical point is obtained by minimizing 
the background held potential. 11 

The inversion symmetry G® implies that there are 


11 The standard proof that the physical point corresponds to the 
absolute minimum of the potential is based on the convexity of 
— In Z( J, A, 11) which is easily shown in the case where the ac¬ 
tion in the functional integral m is real. We note that this is 
not guaranteed in the usual FP quantization procedure since the 
determinant of the FP operator—obtained by integrating out 
the ghost and antighost fields—is not positive definite. We ex- 


V(r,0) = V(—r, 0) = V*(r,0) 
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three more reflexion symmetry axes 12 in the plane (r 3 , rs) 
on top of those identified in Sec. IIII Cl One is the axis 
r 3 = 0 and the two other are the lines which pass through 
the origin and which make an angle ±7 t/3 with the rg = 
0 axis. This is represented in Fig. [T] The origin and 
the points related to it by the translations described in 
Sec. HTTP] now have the symmetry C§ v . 

Using the symmetries of the potential, we can restrict 
the search for its minimum to the elementary cell de¬ 
picted as a shaded triangle on Fig. |T| As discussed in Ap¬ 
pendix [Aj the assumption that charge conjugation sym¬ 
metry is not spontaneously broken at fx = 0 implies that 
the absolute minimum of the background potential lies 
on the axis r$ = 0 in this cell. This is indeed, what 
we find from the detailed investigation of the one-loop 
potential. Furthermore, Eqs. CD and ( III[ i imply that 
£(r, 0 ) = £(r, 0 ) € 1 R on the axis r$ = 0 so that the phys¬ 
ical Polyakov loops, Eqs. Cl and d, are equal and 
real, 


£(H = 0) = £{n = 0) £ R, (72) 

as expected. Indeed, the fact that they are real is con¬ 
sistent with their standard interpretation in terms of the 
free energy of a static quark or antiquark . The fact 
that they are equal simply tells that there is no distinc¬ 
tion between the free energy of a quark and that of an 
antiquark at fj, = 0 , as a result of the charge conjugation 
symmetry. 


B. One-loop results 

To analyze the phase diagram at fi = 0, we track the 
absolute minimum of the one-loop background field po¬ 
tential computed in Sec. HVI as a function of temperature 
for different values of the quark masses. We consider 
the case of two degenerate flavors with mass M u and a 
third flavor with mass M s . In the following, we shall ei¬ 
ther present our results in units of to for dimensionful 
quantities, or compute dimensionless quantities that can 
be compared with lattice results. To get a rough idea 
of scales, a typical value used to fit lattice propagators 



FIG. 2 : The background field potential (in arbitrary units) 
at /i = 0 as a function of r3 for rg = 0 for different tempera¬ 
tures in the degenerate case M u = M s = M. Top figure: M 
slightly larger than the critical value M c ; see Fig. [ 3 ] Temper¬ 
ature increases for curves from bottom to top. The dashed 
line corresponds to the transition temperature. Middle fig¬ 
ure: M = M c with the same conventions as for the top figure. 
Bottom figure: M < M c . The dashed line is for the crossover 
temperature, where the curvature of the potential at the min¬ 
imum is the smallest. 


pect the situation in the present massive extension to be more 
favourable to a convex — In Z (,/, A. ji) since the mass term sup¬ 
presses the large field configurations for which the FP operator 
develops negative eigenvalues. A more thorough investigation of 
these aspects goes beyond the scope of the present work and we 
postpone it to a future work (we stress that this is not particu¬ 
lar to the present model but it is a general issue for background 
field methods). Here, we shall make the conservative assumption 
that in the cases where the background field potential is a real 
function of real variables, the physical point corresponds to an 
absolute minimum. 

12 These are symmetries of the potential, not exactly of the function 
i(r, 0), which gets complex conjugated in the inversion through 
the origin; see Eq. 1711 ) . 


at zero temperature for the SU(3) Yang-Mills theory is 
to « 0.5 GeV [H, [33- 

Depending on the values of the quark masses, we find 
different types of behaviors when we change the temper¬ 
ature, as illustrated on Fig. [D For large masses, the 
absolute minimum presents a finite jump at some tran¬ 
sition temperature, signaling a first-order transition. In¬ 
stead, for small masses, there is always a unique mini¬ 
mum, whose location rapidly changes with temperature 
in some crossover regime. At the common boundary of 
these two mass regions, the system presents a critical be¬ 
havior: there exists a unique minimum of the potential 
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for all temperatures, which however behaves as a power- 
law around some critical temperature. The associated 
nonanalytic behavior of the minimum of the potential as 
a function of temperature is a consequence of the fact 
that, at the critical temperature, the first, second and 
third derivatives in the r% direction vanish, a criterion 
which we used in order to determine the critical line in the 
Columbia plot of Fig. [3] In the degenerate case M u = M s 
we get, for the critical mass, M c /m = 2.867 and, for 
the critical temperature, T c /m = 0.355. We thus have 
M c /T c = 8.07. This dimensionless ratio does not de¬ 
pend on the value of m and can be directly compared to 
lattice results. For instance, the calculation of Ref. [2(ij ] 
yields, for 3 degenerate quarks, (M c /T c ) latt - = 8.32. We 
obtain similar good agreement for different numbers of 
degenerate quark flavors, as summarized in Table |U 


N f 

M c /m 

T c /m 

M c /T c 

(M c /T c ) latt ' 

(M c /T c ) matr - 

1 

2.395 

0.355 

6.74 

7.22(5) 

8.04 

2 

2.695 

0.355 

7.59 

7.91(5) 

8.85 

3 

2.867 

0.355 

8.07 

8.32(5) 

9.33 


TABLE I: Values of the critical quark mass and temperature 
for Nf — 1,2,3 degenerate quark flavors from the present 
one-loop calculation. The values of the dimensionless and pa¬ 
rameter independent ratio M c /T c are compared to the lattice 
results of Ref. [2[| (before to last column). For comparison 
we also show the values from the matrix model of Ref. [24j 
(last column). 

We observe that the critical temperature is essentially 
unaffected by the presence of quarks. It is actually close 
to the one obtained in the present approach for the pure 
gauge SU(3) theory [29}. This is due to the fact that, for 
the typical values of M/T near the critical line, the quark 
contribution to the potential is Boltzmann suppressed as 
compared to that of the gauge sector, as discussed in 
Appendix [B] 

Finally, we compare our results for the ratio M c /T c 
to those of Ref. (24[ based on matrix models. We also 
mention that recent calculations in the Dyson-Schwinger 
approach 0 yield values of the ratio M c /T c that are sys¬ 
tematically smaller than the ones obtained on the lattice. 
The origin of this discrepancy remains to be clarified. 

VI. IMAGINARY CHEMICAL POTENTIAL 

The study of the QCD phase diagram for imaginary 
chemical potential is of great interest. There is no sign 
problem in this case and lattice simulations can be per¬ 
formed in a standard way 00- Thermodynamic po¬ 
tentials at real chemical potential can then be obtained, 
in principle, by analytic continuation. In practice, how¬ 
ever, this is restricted to small values n/T < 1 because 
the theory at imaginary chemical potential has a nontriv¬ 
ial phase structure with nonanalytical behaviors. The 
phase diagram at imaginary chemical potential is also 
interesting in itself as it presents bona fide properties of 



FIG. 3: Columbia plot at n = 0. In the upper right corner of 
the plane ( M U ,M S ), the phase transition is of first-order. In 
the lower left corner, the system presents a crossover. On the 
plain line, the system has a critical behavior. 

QCD 0 , 0 . In this section, we describe the predic¬ 
tions of our one-loop calculation for a purely imaginary 
chemical potential fi = i\Xi- As in the previous case, we 
begin the discussion by analyzing the consequences of the 
symmetry properties discussed in Sec. IIIII 

A. Symmetries 

As before, we seek conditions on the source J and the 
background field A such that the identification A = A, 
with A defined in Eq. CD is consistent. The relations 
CD and CHI imply that the generating function is real 
if we choose both the source and background field to be 
real. We have, in that case, 

Z(J,A,iHi) = Z(-J,-A , -inf) = Z*(J,A,ifj,i). (73) 

The average field A of Eq. CD is thus real and can be 
consistently identified with A. The background field po¬ 
tential V(r, inf is also real, the Polyakov loop functions 
are related by complex conjugation, £(r, inf) = £*(r, inf), 
and so are the physical Polyakov loops: 

(74) 

The Polyakov loop is complex for generic values of j±i and 
its argument actually plays the role of an order parameter 
for the various transitions described below [46[. Finally, 
the background field potential being a real function of 
real variables as in the case /r = 0, we adopt the same 
prescription as before to choose the background field for 
the evaluation of physical observables, namely we mini¬ 
mize the potential V. 

The symmetry properties of the potential V(r,inf for 
a generic value of p, are summarized in Fig.d The white 
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FIG. 4: The white (black) dots are equivalent up to sym¬ 
metry transformations that exist for a generic chemical po¬ 
tential. In particular, the system is invariant under reflexion 
with respect to the plain lines. The shaded equilateral trian¬ 
gle corresponds to the elementary cell that is repeated all over 
the plane. For imaginary chemical potential /z = i/Zi, there 
exists specific values, Hi = mv//3, with n £ Z, for which the 
theory is invariant under inversion about the origin and the 
white and black dots are equivalent. The shaded triangle is 
symmetric under reflexion about the dashed vertical line so 
that the elementary cell is now half of the shaded equilateral 
triangle. Finally, the symmetries of the effective potential im¬ 
ply that the elementary triangle of the theory at imaginary 
chemical potential Hi ±27t/(3 /3) can be obtained from that of 
the theory at Hi by a rotation of ±27r/3 around its center. 


dots are all images of each others by the global color sym¬ 
metries m and the gauge transformations (l42l) and are 
thus all physically equivalent. The same holds for the 
black dots but, contrarily to the case Hi = 0 discussed 
in the previous section, the black and white dots are not 
equivalent. The elementary cell, to which one can restrict 
the analysis, is thus twice as large as that of the previous 
section. It can be chosen as the shaded equilateral tri¬ 
angle shown in Fig. Q] Observe that the origin [and the 
equivalent points obtained by the translations (l42l) ] have 
the symmetry C^ v . 

Remarkably, the elementary cell has an extra symme¬ 
try when the imaginary chemical potential Hi is a mul¬ 
tiple of 7 t/(3 j3). Using the general periodicity property 
V (r, /z) = V(r, fj, + 2in/ft), one can reduce the discussion 
to the interval 0 < Hi < 2n//3. The case /z, = 0 has been 
discussed before and possesses a Cq v symmetry around 
the origin [and all the equivalent points obtained by the 
translations (S3)]. This is also true for Hi = n /j3. Indeed, 


using Eq. SH) and the 27 r// 3 -periodicity in /z,, we have 

V(r,in / 0) = V(—r, —in/ (3) = V(—r,in/ft). ( 75 ) 

This corresponds to an inversion about the origin and 
we thus retrieve the same set of symmetries as in the 
case in = 0, depicted in Fig. [TJ Finally, the cases 
Hi = 2nir/(3ft) and Hi = (2n + l)7r/(3/3), with n integer, 
can be obtained from the cases Hi = 0 and Hi = ft, 
respectively, by using the relation (|46l) . The symmetries 
are similar to those depicted in Fig. [T| where, however, 
the point with Cq v symmetry is not the origin anymore 
but one of the two other vertices of the elementary tri¬ 
angle [and their images by the translations (l42l) ]. For 
Hi = 7t/(3 ft) and Hi = 47T/ (3/3), this is the upper ver¬ 
tex, located at (2tt, 27t/>/ 3), whereas for Hi = 27r/(3/3) 
and Hi = 57t / (3/3), it is the lower vertex, located at 
(27t, — 2-7r/v / 3). Using the axes of symmetry correspond¬ 
ing to global color transformations (l3ijl) (the plain lines in 
Fig. CD, one sees that the elementary triangle is rotated 
by ±27t/3 each time the imaginary chemical potential is 
shifted by ±27r/(3/3). The corresponding Polyakov loop 
acquires a phase ±27t/3. 


B. One-loop results 

We now present our results for the phase diagram 
in the plane {Hi/T, T). Using the property (ISTl) . it is 
sufficient to study the interval 0 < Hi/T < 27r/3 and 
the phase diagram is symmetric around Hi/T = 7r/3. 
For simplicity, we consider the case of three degenerate 
quarks with M u = M s = M. For any mass larger than 
M c (h = 0) — 2.8 m, such that the transition at vanishing 
chemical potential is first-order (see Fig. 0 ), we find that 
the transition persists at nonvanishing Hi, as depicted in 
the top panel of Fig. El Two first-order phase transition 
lines join at the symmetry axis Hi/T = 7r/3. Moreover, 
for sufficiently large temperatures, there exists a first- 
order phase transition across Hi/T = 7r/3, where the ar¬ 
gument of the Polyakov loop has a finite jump, signaling 
the spontaneous breaking of the Roberge-Weiss symme¬ 
try 13 [46[. In the same range of temperatures, we also 
observe that the reflection symmetry of the potential at 
Hi/T = 7r/3 with respect to the median passing through 
the vertex ( 27 t, 27 t/v / 3 ) is spontaneously broken. As Hi/T 
changes from 7 t/ 3 _ to 7r/3 + , the minimum has a finite 
jump from one to the other side of this median. This 
leads to a cusp in the partition function at Hi/T = 7r/3, 
as anticipated in [it|- The three first-order lines meet 
at a triple point, where the system can coexist in three 
different phases characterized by different values of the 
argument of the Polyakov loop or, equivalently at one- 


13 When the Roberge-Weiss symmetry is not broken, Eq. (1521) im¬ 
plies that e~ l7r / 3 £(i7rT/3) G R so that Arg^(z 7 rT/ 3 ) = 7 t/3 . 
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FIG. 5: Phase diagram in the plane (/ u/T,T) for different 
values of the degenerate quark mass M. Top panel: M = 3 m 
is larger that the critical mass at vanishing chemical potential 
M c ( 0) ~ 2.8 m. Middle panel: M = 2.6 m is below M c ( 0) but 
larger than the second critical mass M c (iirT/3) ~ 2.2m; see 
text. Bottom panel: M = 2m is smaller that M c (inT/3). 
Plain lines correspond to first-order phase transitions, dashed 
lines to crossovers, black dots represent Z 2 critical points, and 
the empty circles are triple points. 


loop order, by the value of the background field in the 
rs-direction at the minimum of the potential. 

For M = M c (n = 0), the transition at vanishing chem¬ 
ical potential is second order and there appears, in the 
(fj,i/T,T) phase diagram, a couple of Z 2 critical points 
which terminate the lines of first-order transitions de¬ 
scribed above at m/T = 0 and m/T = 2n/3. Decreas¬ 
ing the mass M further, these critical points penetrate 
deeper in the phase diagram towards m/T = 7t/ 3, as 


Arg i 



FIG. 6: Argument of the Polyakov loop in the (m/T, T) plane 
for a degenerate quark mass M in the range M c (inT/3) < 
M < M c (0), corresponding to the middle panel of Fig. O first- 
order transitions appear as a discontinuities in this quantity. 
The discontinuity vanish at the Z 2 critical points. At the 
triple point on the line m/T = 7r/3 , the argument of the 
Polyakov loop can take three different values. 


shown in the middle panel of Fig. 0 At a critical value 
M c (n = inT/3) ~ 2.2m, the two Z 2 critical points merge 
at the symmetric point m/T = 7t/ 3 and give rise to a 
tricritical point which terminates the vertical line of first- 
order transition [49j. The horizontal lines of first-order 
transitions for m/T ^ 7t/ 3 have completely disappeared 
and are replaced by crossovers. For M < M c (*7tT/ 3), the 
picture is the same with, however, the tricritical point 
ending the first-order transition line at m/T = 7r/3 re¬ 
placed by a Z 2 critical point, as shown in the lower panel 
of Fig. [5l As a further illustration, we display in Fig. [G] 
the argument of the Polyakov loop in the (m/T, T) plane 
for the intermediate mass case, M c (i7r/3) < M < M c ( 0). 
Similar plots can be made in the other cases as well. We 
mention that similar results were obtained using effective 
matrix models in Ref. . 

The tricritical point at m/T = 7t/ 3 and M = 
M c (inT/3) features particular scaling laws which con¬ 
strain the ^q-dependence of the critical quark mass and 
temperature, M c and T c , at the Z 2 critical points. If 
the scaling window is broad enough, this may put con¬ 
straints on the phase diagram away from the tricritical 
point and possibly for real chemical potential as well [49j . 
Our results can be qualitatively described by the follow¬ 
ing model potential 

- H (76) 

where a oc M c (inT/3) — M, b is a decreasing func¬ 
tion of the temperature, and the ^-breaking term h oc 
(7t/3) 2 — (m/T) 2 characterizes the explicit breaking of 
the Roberge-Weiss symmetry, that is, the distance to the 
symmetry axis m/T = 7t/3. Here, <j> parametrizes the 
curve in the (rz,rs) plane along which the minimum of 
the potential V(r, i\x/) moves. Equation (1761) is the mini¬ 
mal model for describing the tricritical scaling. For h = 0 
and a < 0, it describes a first-order transition with a 
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FIG. 7: Top panel: The quark mass at the critical point as a 
function of x = (7r/3) 2 — (/ii/T) 2 . The white squares are the 
results of the present one-loop calculation and the plain line 
is the fit 6.15 + 1.85a: 2 ' 5 . Bottom panel: The temperature in 
units of m at the critical point as a function of x. The plain 
line is the fit 0.35 + 0.012 x 2/5 - 0.0055 a: 4/5 . 


placing /ii —> —ifj,. In particular, they can be compared 
to direct calculations of the phase diagram at real chem¬ 
ical potential, as discussed in the next section. 

It is also interesting to compare our results to lattice 
calculations. Following Ref. [49(, we write 


M c (n) 

W 



(77) 


A fit of our results at /i = yields M tr i c ./T tr i c . = 6.15 
and K = 1.85, to be compared with the lattice result 
of Ref. 0, (Mt ric ./Tt r i c .) latt - = 6.66 and X latt - = 1.55 
for 3 degenerate quark flavors. 14 As for the case /r = 0, 
we note that the values of the fitting parameters obtained 
in the recent Dyson-Schwinger calculation of Ref. [HI are 
systematically smaller than the one obtained here and on 
the lattice. Again, the origin of this discrepancy remains 
to be elucidated. 


VII. REAL CHEMICAL POTENTIAL 

The case of real chemical potential is where the lat¬ 
tice methods suffer from a severe sign problem because 
the QCD action in the functional integral measure is not 
real. As we shall discuss below, this also affects con¬ 
tinuum approaches, although in a much less severe way. 
As before, we first discuss the symmetries of the prob¬ 
lem and the issue of choosing the source and background 
field in a consistent way. This is where the remnant of 
the sign problem appears. We then present the results of 
our one-loop calculation. 


jump in 4> at the minimum oc y/—a. Instead, the phase 
transition is continuous with a Z 2 critical point for for 
a > 0 and a tricritical point in the limiting case a = 0. 
This qualitatively describes the various cases depicted in 
Fig. [5]on the axis ^i/T = 7t/3. 

Now, for h 0, the potential (1761) predicts a crossover 
as a function of temperature for a > 0 and either a 
first-order transition for sufficiently large negative a or 
a crossover for small negative a. This corresponds to the 
various cases of Fig. [5] away from the axis Hi/T = 7r/3. 
The critical point in the middle panel is the limiting 
case between the first-order transition and the crossover 
for a < 0. The corresponding critical values scale as 
a c ~ h 2 / 5 and b c ~ / 1 4 / 5 . This directly translates to the 
scaling law for the critical quark mass and the critical 
temperature. Such scalings are well reproduced by our 
results, as shown in Fig. [3 The scaling of the critical 
quark mass directly follows that of the parameter a in 
the above model. However, the parameter b is not only 
a function of the temperature but also depends on the 
mass. As a consequence, the scaling of the critical tem¬ 
perature has both a contribution in ft 2 / 5 and a contribu¬ 
tion in / 1 4 / 5 . As mentioned before, such scaling relations 
can be extended to arbitrary chemical potential by re- 


A. Symmetries 

The main difficulty that we have to face at real chemi¬ 
cal potential is the fact that the generating functional J5]) 
is not real in general; see Eq. C51) . As a consequence, the 
average value A of the gluon field, obtained from Eq. (fTIT) . 
is typically not real and one needs to consider complex 
background fields 15 (and possibly complex sources) in or¬ 
der to be able to identify A = A in a consistent way. How¬ 
ever, one has then to deal with a potential V (r, /i) which 


14 We mention that the value A' latti = 1.55 of Ref. I20II corresponds 
to the case Nf = 1. Although these authors do not mention 
the corresponding value for Nf =3, we observe that keeping 
the same value gives a reasonable estimate of the ratio M c /T c at 
/r = 0 quoted in that reference. 

15 With a complex background A the meaning of the shift © un¬ 
der the path integral needs to be reconsidered. In fact there 
is no contradiction because, strictly speaking, in the derivation 
of the loop expansion formula, the shift needs to be done after 
the original (real) contours of integration for the components of 
the gauge field have been deformed to pass through the relevant 
saddle point. At this level, the variables of integration are a pri¬ 
ori complex and 0 is a mere complex change of variables for 
contour integrals. 
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is a complex function of complex variables. Not only is 
the problem more intricate but it is not even clear what 
is the correct procedure to identify the physical value 
of the background field needed to compute observables 
(which should end up being real despite the fact that the 
extremum of V ( r , /x) is possibly complex). For instance, 
these do not correspond in an obvious way to the absolute 
minimum of the potential anymore and might instead be 
saddle points. In that case, it is not clear which saddle 
point is the preferred one. 

The problem is greatly simplified—although not com¬ 
pletely solved—by choosing the source and background 
field as 

j =(5)' M^ 8 )- (78) 

with J 3,8 and A 3,8 real, such that 16 

RiJ = -J * and R X A = -A*, (79) 

where R± is the reflexion with respect to the 8-axis, de¬ 
fined in Eq. (1511) . Using Eqs. (1151) and (1551) . one sees that 
the generating function is real, 17 

Z[J, A , fj] = Z(RiJ, Ri_A, g) = Z*(J, A, a«). (80) 

This has two important consequences. First, the aver¬ 
aged gluon field (fTTl) is of the form A = (A 3 ,*A 8 ) with 
A 3,8 real, which guarantees that the identification A = A 
is consistent. Second, it implies that the effective poten¬ 
tial is a real function of the real variables r 3 , 3 = fig A 3,8 : 

V(r, At) = V (r 3 , ir 3 , aO = V*(r 3 ,ir 8 , g). (81) 

Note that the potential is not invariant under transla¬ 
tions and reflexions in the (imaginary) r 3 direction but 
these symmetries are still valid in the r 3 direction. Con¬ 
sequently, it is sufficient to study the band r3 S [0, 2-k\. 

Now comes the question of selecting the correct ex¬ 
tremum. We first observe that in the case g = 0, both 
U(r3,r8,0) and V(r 3 ,i r s,^) are real for real r 3 ^- In the 
former case, the relevant extrema are the absolute min¬ 
ima of the potential, which lie on the r 3 = 0 axis, as 
discussed in Sec. IV Al It is clear that if we now change 
r$ —> ir 3 , the curvature in the r% direction changes sign 
and these minima turn into saddle points. For g ^ 0 not 


16 A similar proposal has been made in Ref. ;27( . where the authors 
also propose to search for saddle points of the effective potential 
V(r 3 , irs, //) in the plane ( 7 * 3 , rs); see below. 

1 ‘ By using the color transformations of Sec. IIII Cl one can find 
other ways of having a real generating function where both r 3 
and rg are complex. These are however equivalent to the choice 
considered in the main text and lead to identical physical pre¬ 
dictions. We have also investigated the case where both r 3 and 
rg are purely imaginary, for which V(j\ ji) is also real. However, 
the one-loop expression of the potential does not make sense in 
that case. 






h 

FIG. 8: Contour plots of the background field potential 
V(r 3 ,ir 8 , g) in the plane (r 3 ,rg) for M u = M s = 4 m and 
g = 0.6m for T < T c (top) and T > T c (bottom), with 
T c = 0.36m. The saddle point where the potential is the low¬ 
est is indicated by the black dot. We clearly see the jump 
at T c and the fact that the saddle points sit at rs A 0. The 
dashed curve is an approximate analytic expression of the po¬ 
sition of the saddle point in the rs-direction as a function of 
its r 3 coordinate; see Appendix m 


too large, we can follow these saddle points, which con¬ 
tinuously move away from the r% = 0 axis. Which one to 
choose as the physical point is, however, not clear. Un¬ 
like in the case of a real action, we have no first principle 
argument for preferring a particular extremum. This can 
be seen as a remnant of the sign problem of lattice sim¬ 
ulations; see (50| f° r a similar discussion in the context 
of a charged scalar field. Even if the complex nature of 
the action is not an obstacle to perform calculations using 
continuum approaches, it leads to ambiguities in selecting 
the physical solution out of the various extrema. Here, 
we always choose (arbitrarily) the saddle point for which 
the real function V(r 3 ,i r s, t) is the lowest. Although 
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FIG. 9: Columbia plot for real chemical potential. The first- 
order region contracts with increasing /r. The dotted line 
corresponds to fi = 0. Successive plain lines correspond to a 
chemical potential increased by steps Sfx = 0.1 GeV. 

this seems reasonable at small /r by continuity with the 
case fi = 0, we have no guarantee that this is a valid 
procedure—if at all -for all values of fi. 

Let us finally mention that with the choice of back¬ 
ground field coordinates considered here, the Polyakov 
loops are always real, as follows from the relations (l30l) 
and (1371) : 

Z(r, n) = £(Rir, n) = t*{r, n) ( 82 ) 

and similarly for £(r,n), where r = (r^,ir^). As for the 
case fi = 0, the fact that the Polyakov loops are real 
is consistent with their interpretation in terms of the 
free energy of static quarks or antiquarks. In general, 
t{r,n) ^ £{t,h) at fi yf 0, which simply reflects the fact 
that the free energy of a static quark differ from that of 
its antiparticle at finite ji. We note though that at tree- 
level we have the additional relation £(r,fj) = £(—r,fi) 
[see Eq. flSZD] so that the two functions coincide on the 
axis rs = 0. There is, however, neither any reason of 
symmetry nor physical argument for this to be the case 
and we do not expect it to be true at higher orders. 


B. One-loop results 

We concentrate here on the determination of the crit¬ 
ical line in the (M U ,M S ) plane. We perform a similar 
analysis as in Sec. El but this time identifying saddle 
points. A typical situation is illustrated in Fig. [51 where 
we see the jump of the deepest saddle point at the first- 
order transition. We recover the fi = 0 results and we 
find that the region of first-order transition shrinks to¬ 
wards large quark masses when the chemical potential 
increases, as shown in Fig. [9] 



FIG. 10: The critical mass M c (ix)/T c (n) as a function of x = 
(n /3) 2 + (fi/T) 2 for the case of degenerate quark masses M u = 
M s . The crosses are the result of the direct calculation. The 
line is obtained from the tricritical scaling law ([13; see also 

Fig. m 


We plot the /x dependence of the critical quark mass 
in the degenerate case M u = M s in Fig. [TO] and we com¬ 
pare it with the expectation from the tricritical scaling 
Gil) extrapolated from the imaginary chemical potential 
region. We see that this describes well the data at real 
chemical potential. 

In our analysis, we observe that the location of the sad¬ 
dle point is typically at r§ ^ 0, which induces a difference 
between the averages of the Hermitic conjugate Polyakov 
loops £{\x) and £(fx) already at tree-level. As already em¬ 
phasized, this is physically related to the different free 
energies of static quarks and antiquarks at finite fi. An 
example is depicted in Fig. 1111 which shows the temper¬ 
ature dependence of the averaged Polyakov loops in the 
region of first-order transition. We observe a significant 
difference between £(fx) and £(fx) below the transition 
temperature, whereas they essentially agree in the high- 
temperature phase. In other words, the energetic price 
to pay for a static antiquark is much higher than that for 
a quark (at fi > 0) in the quasi-confined, low tempera¬ 
ture phase, where the Polyakov loops are small, whereas 
it is essentially the same in the high-temperature decon- 
fined phase. This can be understood analytically from 
an approximate calculation presented in Appendix iBl 

We conclude this section by mentioning that previous 
studies of the phase diagram with background field meth¬ 
ods in the functional renormalization group and Dyson- 
Schwinger approaches of Refs. [12, EH have employed an¬ 
other criterion than the one used here to determine the 
physical properties of the system. Instead of searching 
for a saddle point of the function V (7*3, irg, /j.) as we pro¬ 
pose here, the authors of Refs. [HI, [l(| define the physical 
point as the absolute minimum of the function V (r3, 0, n) 
as a function of . We have repeated our analysis using 
this procedure for comparison. A clear artifact of this 
procedure is that on the axis r% = 0, the tree-level ex¬ 
pressions of the Polyakov loops £(/i) and £([j) are equal, 
as already mentioned. However, we have found that both 
criteria give essentially the same critical temperatures in 




16 



T/m 



FIG. 11: Upper panel: Average of the traced Polyakov loop 
and of its Hermitic conjugate as a function of temperature 
for degenerate quark masses M u = M s = 4m and chemical 
potential p = 0.6m, as in Fig. [8] At low temperatures, the 
average of the Polyakov loop (higher curve) differs from that 
of its Hermitic conjugate (lower curve). The dotted curve 
corresponds to the determination of the order parameter cal¬ 
culated for V8 = 0. The lower panel shows the static quark 
and antiquark free energies obtained from the Polyakov loops 
of the upper figure as i = exp(—/3F q ) and £ = exp(— /3F q ). 


our calculation. This is illustrated in Fig. [TTJ 

That the critical temperatures are not significantly 
modified in these two prescriptions can be traced back 
to the relative smallness of the values of r§ obtained by 
following the saddle points in our procedure; see Fig. [8] 
This, in turn, originates from the strong Boltzmann sup¬ 
pression of the (heavy) quark contribution to the poten¬ 
tial, which is responsible for the departure of the saddle 
point from the axis r§ = 0, as discussed in Appendix iBl 
We point out that the situation mig ht be very differ¬ 
ent in the case of light quarks [12, lL 6| and that different 
procedures for identifying the relevant extremum of the 
potential may have more dramatic consequences. This 
needs to be investigated further. 


VIII. CONCLUSIONS 

We have studied the phase diagram of QCD with heavy 
quarks at leading perturbative order in a massive exten¬ 


sion of the Landau-DeWitt gauge in the context of back¬ 
ground field methods. The richness of the phase diagram 
is reproduced for both real and imaginary chemical po¬ 
tential and we obtain parameter free values for dimen¬ 
sionless ratios of quark masses over temperature at crit¬ 
icality which agree well with results from lattice simula¬ 
tions. This extends previous studies of the confinement- 
deconfinement phase transition in pure Yang-Mills theo¬ 
ries 013 and adds to the list of lattice results that 
can be—at least qualitatively if not quantitatively— 
described bv this modified (massive) perturbative scheme 

[nafflsi. 

A natural extension of the present work is to com¬ 
pute the next-to-leading perturbative corrections both 
to the potential and to the order parameter and to study 
whether the results converge toward lattice values. This 
can be done along the lines of Refs. [30j, [U]. Finally, it 
would be of interest to investigate the phase diagram in 
the light quark region. However, this requires the treat¬ 
ment of chiral symmetry breaking and involves some re¬ 
finement of the present approach. We defer such studies 
for future work. 
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Appendix A: Charge conjugation symmetry at p = 0 

In the presence of a nonvanishing background field, one 
must be careful in assessing whether a symmetry is (ex¬ 
plicitly or spontaneously) broken since the background 
field itself explicitly breaks some symmetries, e.g., the 
global color or the charge conjugation symmetries. To 
probe a possible symmetry breaking, one must consider 
a gauge-invariant observable, which is independent of the 
background field. In the present approach, any such ob¬ 
servable O can be obtained as the value of some function 
at the absolute minimum 18 r m i n (/x) of the poten¬ 
tial: 

O(fJ-) = 0(r mi n (/j),/j). (Al) 

Let us for instance consider the case of the charge con¬ 
jugation symmetry C , discussed in the main text. A C- 
odd observable shares the same symmetry properties as 
the potential V(r,/i) (such that, again, one can restrict 
to the elementary cell in the Cartan plane) except that 

0(r, M ) = -0(-r,-/z). (A2) 


18 We consider the case where the physical point is the absolute 
minimum of the potential, relevant for the cases = 0 or /iG iR 
discussed here. 
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In particular, for /i = 0, we can use the color symmetry 

(ED to get 

O(r 3 ,r s ,0) = -0(r 3 ,-r 8 ,O), (A3) 

which implies that 0 (r, 0) = 0 on the axis r 8 = 0 and, 
by symmetry, on any of the dashed axes in Fig. [Q Now, 
assuming that the charge conjugation symmetry is not 
broken at fi = 0, we have 0(0) = 0(r m i n (O), 0) = 0 
for any such observable. In other words, the functions 
0 (r, 0) which are odd with respect to the mirror symme¬ 
tries mentioned above (i.e. r 8 —> —rg in the elementary 
cell) must vanish at the minimum of the potential. For 
this to be true for all odd functions of r 8 (we assume that 
any such function can be reached from the operators of 
the theory), the extremum must necessary lie on the axis 
rg = 0. This shows that 

C manifest => r 8 = 0 r 8 ^0=>C broken. 

(A4) 

The above conclusion can be strengthened if we make 
the (reasonable) assumption that every minimum of the 
potential V (r, y) in the elementary cell of Fig. [4] corre¬ 
sponds to a distinct physical state. In that case we con¬ 
clude that, for any value of /z, 

C manifest •£=> r 8 = 0. (A5) 

In particular this implies that the relevant extremum of 
the potential must lie at rg ^ 0 at y 0 since C is 
explicitly broken. 

Appendix B: Large mass approximate formula 


with the tree-level Pokyakov loops in the fundamental 
representation 

£f(t) = ^ 22 eirp = 7> + 2e*27S cos -y ^ (B5) 

p 

and £f(t) = £f(~ r). 

In the range of temperatures and quark masses con¬ 
sidered in this study, the Boltzmann suppression factors 
are /(/3m « 3) ~ 5.6 x 10~ 2 for the massive gluon modes 
and f((3M > 6) < 6.2 x 10 -3 for quarks. Note that 
in the limit of very heavy (nonrelativistic) quarks, one 
has f((3M) ~ (/3M) 3 / 2 e“^ M /(-\/27r 3//2 ). This gives a rea¬ 
sonably good estimate for the cases considered here, with 
/3 M > 6. This relative suppression of the quark contribu¬ 
tion as compared to that of gluons explains the fact that 
the critical temperatures obtained in the present work 
are essentially insensitive to the presence of quarks; see 
Table U 

The relative suppression of the quark contribution also 
explains the smallness of r 8 at the extremum of the po¬ 
tential V(rg,ir s , /z) obtained in Sec. IVIII Indeed, as dis¬ 
cussed in Sec. El the extremum of the potential always 
lies on the axis r 8 = 0 at /z = 0. A nonzero value of r 8 
at /z y^ 0 is, thus, entirely due to the quark contribution 
and is thus controlled by the suppression factor f((3M). 

It is an easy exercise to compute the actual value of r 8 
as a function of rg at the extremum at leading order in 
f((3M). Writing 


For the range of temperatures and quark masses stud¬ 
ied here, we have typically (3M > 6 and (3m ss 2.75. 
The massive degrees of freedom are thus essentially de¬ 
scribed by classical Boltzmann statistics [H, [28j . Using 
exp(— (3e q ) -C 1 in the integral (1CT1) . the pure gauge po¬ 
tential (l59l) can be approximated as 


1^5# « - 12/G8m)4,(r), (Bl) 

where we defined the function 

1 r°° - t 2 

/Or) = — / dyy 2 e-Vv^ = — K 2 {x), (B2) 

7T 2 J 0 7T 2 

with K v (x) the modified Bessel function of the second 
kind and where 


£a(t) = - '22 e= - (1 + cos rg + cosr + + cos r_) 


(B3) 

is the tree-level Polyakov loop in the adjoint representa¬ 
tion, with r± = (rg ± rg\/3)/2. 

Similarly, the quark contribution (1571) is approximated 
as (we consider Nf degenerate flavors with mass M) 


« -3 N f f(fiM) [e-^Mr) + e^t F (r)\ (B4) 


V(r 3 ,ir 8 ,n) = V(r 3 , ir s , 0) + SV(r 3 , ir s , /z), (B6) 


and using the fact that d rg V(r 3 , irg , 0)| rs= o = 0 V r3, one 
shows that, at leading order in SV, the value of r 8 at the 
extremum is given by 


r e g xt -(r 3 ) 


d r8 dV(rg,ir 8 ,n) 

dr a V(r 3 Jr s ,0) 


(B7) 


Using the approximate expressions for the potential ob¬ 
tained above, a simple calculation yields 


r! xt -(r 3 ) 


16Nf 


f((3M) sinh(/3/z) sin 2 (r 3 /4) 

1 “ “ n ) 2 + 18 /(/3 m ) cos(r 3 /2) 


(B8) 

We have plotted this formula against the r 3 coordinate of 
the saddle point in Fig. [HI which shows that it describes 
well the results from the complete calculation. Inciden¬ 
tally, Eq. (1B8I) also explains the relative smallness of the 
ratio (£ — £)/(£ + £) in the high-temperature phase as 
compared to the low-temperature phase in Fig. 1111 This 
originates from the rapid increase of r| xt ’(r 3 ) with in¬ 
creasing r 3 in the range of interest; see Fig. [S] 
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